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The  Assimilation  of  Marine  Surface  Data  into  the  Navy 
Operational  Global  Atmospheric  Prediction  System 


ABSTRACT 

numerical  models  of  the  ocean  and  the 
atmosphere  become  more  sophisticated,  it 
becomes  necessary  to  account  for  the  inter¬ 
change  of  energy  at  the  air/sea  interface. 
Meteorological  modeLs  require  accurate  esti¬ 
mates  of  sea  surface  temperature  and  ice 
fields.  Oceanographic  models  require  accurate 
estimates  of  the  surface  winds  and  heal  flu.res. 
Thus,  it  is  beneficial  for  scientists  in  each  dis¬ 
cipline  to  have  some  level  of  understanding 
and  appreciation  for  how  these  various  param¬ 
eters  are  obtained. 

With  that  in  mind,  we  tvill  describe 
how  a  meteorological  analysis/forecast  system 
uses  rnaritte  data  to  estimate  the  surface  rcirtd 
field,  icith  particular  emphasis  on  the  arralysis 
proce.ss  and  the  resultirrg  nowcast  product.  Tire 
Navy's  data  analysis  system  uses  the  mrrlti- 
var-iate  optimum  interpolation  technique  to 
produce  updates  to  backgrourrd  fields  provided 
by  the  forecast  model.  The  advarrtage  of  this 
technique  is  that  the  relative  weights  giverr  the 
background  and  the  obsenrations  are  deter¬ 
mined  by  taking  into  account  the  estimated 
error  properties  of  each.  i4s  rretv  data  sources 
become  available,  the  appropriate  statistics 
defining  the  errors  of  the  observations  must  be 
determined.  Since  we  have  recerrtly  added  .sev¬ 
eral  'nerv  data  sources  to  the  a.ssirnilatiorr  cycle, 
new  values  of  the  observation  errors  harte  been 
determined  for  all  marine  surface  data.  The 
improvements  made  by  these  nerv  data  source.s, 
ittcluditrg  special  serrsor  microivave/irrrager 
tvind  speeds  and  syrrthetic  tropical  obserro- 
tions,  rvill  alsit  be  illustrated. 

INTRODUCTION 

iven  the  current  emphasis  on  global  climate 
change,  the  scientific  community  has  begun 
to  look  more  closely  at  the  total  global  system 
iuid  the  relationships  among  its  various  compo¬ 
nents.  As  a  result,  research  interest  in  coupled 
gkibal  oo'iui/atmosphere  motlels  has  increased 
dramatically.  As  we  study  and  define  the  inter- 
<  liangi'  of  energy  at  the  air/oc^aT’  interface,  sci¬ 
entists  from  each  dist  itiiine  must  also  develop  a 
better  understanding  of  how  the  parameters  of 
the  comiilementary  system  are  derived. 

Although  the  Navy  is  not  directly  involved  in  the 
climate  (>ffort,  they  have  long  been  at  the  fore¬ 
front  of  air/sea  research  to  support  the  re(|tiire- 
rnents  of  their  operations  in  the  marine  bound¬ 
ary  layi'r.  These  efforts  include  but  are  not  limit¬ 
ed  to  numerical  modeling  of  the  atmosphere  and 
ocean 


Accurate  nowcasts  and  forecasts  of 
atmospheric  conditions  are  essential  to  the  Navy 
for  several  reasons.  Most  obviously,  the  atmo¬ 
sphere  affects  ship  and  aircraft  routing,  aircraft 
carrier  operations,  and  weapons  systems  perfor¬ 
mance.  Not  to  be  overlooked,  however,  is  the 
impact  of  tiie  atmosphere  on  the  oceanographic 
products  produced  by  the  Navy.  Currently, 
atmospheric  winds  and  heat  fluxes  provide  forc¬ 
ing  for  ocean  thermodynamic,  hydrodynamic, 
wave,  and  ice  models  (Clancy,  1987),  while  sea- 
surface  temperature  and  ice  analyses  improve 
the  boundary  layer  parameterizations  within  the 
meteorological  forecast  model.  Through  this 
interactive  coupling,  even  the  traditional 
oceanographic  domains  of  search  and  rescue, 
acoustics,  and  arctic  submarine  operations  are 
directly  linked  to  the  atmosphere. 

Various  meteorological  and  oceano¬ 
graphic  products,  such  as  those  mentioned 
above,  are  produced  and  disseminated  by  the 
Navy’s  Fleet  Numerical  Oceanography  Center 
(F’NOC),  while  key  research  and  development 
tasks  are  performed  by  divisions  of  the  Naval 
Research  Laboratory  (NRL)  in  Monterey, 
California,  co-located  with  FNOC,  and  at  Stennis 
Space  Center  in  Mississippi.  While  achievements 
in  air/ocean  model  coupling  are  driven  by  Fleet 
requirements,  these  efforts  are  supported  by  the 
most  contplete  global  database  in  the  world.  In 
addition  to  conventional  radiosonde  and  surface 
data,  FNOC’s  meteorological  database  includes 
satellite-derived  temperature  soundings  and 
satellite-derived  wind  information,  along  with 
synthetic  observations  in  the  vicinity  of  extra- 
tropical  and  tropical  storms.  Although  much  of 
this  data  is  available  over  the  Global  Tele- 
conununications  System  (GTS),  P'NOC  was  the 
first  operational  center  to  have  real-time  access 
to  the  marine  surface  wind  data  from  the  special 
.s('nsor  microwave/imager  (SSM/1)  flown 
onboard  the  Defense  Meteorological  Satellite 
Program  spacecraft.  While  the  SSM/1  data  have 
recently  Ireen  made  available  to  the  National 
Meteorological  Center,  their  distribution  and 
operational  ttse  are  still  limited  at  this  time. 
Furthermore,  classified  observations  and  locally 
generated  synthetic  observations  are  not  avail¬ 
able  to  outside  users. 

To  take  full  advjmtage  of  this  database 
for  improved  meteorological  analysis  (nowcast ) 
anil  forecast  products,  significant  upgrades  wen' 
made  to  the  Navy  Operational  Global  Atmo¬ 
spheric  Prediction  System  (NOGAI’S)  in  .lanuarj 
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of  1988.  Like  most  global  data  assimilation  sys¬ 
tems,  NOGAJ^S  consists  of  four  components  that 
are  inextricably  connected — data  quality  con¬ 
trol,  daUi  analysis,  model  initialization,  and 
model  forectast.  The  newer  generation  of 
NOGAPS  included  mtyor  upgrades  to  all  of  these 
c'omponents.  ITie  grid-point  forecast  model  was 
replaced  with  a  global  spectral  model  (Hogan 
and  Rosmond,  1991),  the  successive  corrections 
analysis  sclieme  was  replaced  by  a  three- 
dimensional  multivariate  optimum  interpolation 
(MVOI)  analysis  (Barker,  1992;  Goerss  and 
Phoebus,  1992a),  the  balance  equation  was 
replaced  by  a  nonlineai  iioinial-mode  initializa¬ 
tion  (Hogan  et  al.,  1991),  iuid  significant 
improvements  were  made  to  the  data  quality 
control  (Baker,  1992).  Since  1988,  modifications 
such  as  higher  resolution,  better  physical  param- 
eterizations,  and  improved  data  handling  have 
continued  to  increase  the  skill  of  the  analysis 
;md  forecast  products.  As  a  result,  NOGAPS  is 
now  competitive  with  the  global  data  assimila¬ 
tion  systems  run  at  other  operational  centers. 

This  paper  will  concentrate  on  the 
near-surface  analyses  over  the  oceans — an  area 
of  obvious  interest  to  the  Navy,  and  an  area 
where  F'NOC’s  unique  database  can  be  fully 
exploited.  The  following  sections  describe  this 
database  and  the  MVOI  analysis  technique  in 
more  detail.  The  final  section  will  also  discuss 
the  impact  of  new  data  sources  on  the  analyzed 
products.  We  hope  to  provide  the  oceanogra¬ 
pher  with  an  overview  of  the  process  that  takes 
meteorological  data  from  a  number  of  sources 
that  are  randomly  distributed  in  space  and  time 
and  from  them  produces  spatially  coherent  grid- 
ded  wind  fields  in  the  marine  bouitdaiy  layer. 

DATA  ANALYSIS 

he  global  analysis  product  is  produced  for 
sixteen  pressure  levels  (the  standard  levels 
from  1000  millibars  (mb)  to  10  mb)  on  the 
Gaussian  grid  of  the  NOGAPS  spectral  forecast 
model.  The  Gaussi2m  grid  points  are  separated 
by  approximately  1.5',  but  are  thinned  near  the 
poles  for  the  data  analysis.  The  analyzed  meteo¬ 
rological  variables  are  geopotential  height  k  and 
the  zonal  aitd  meridional  wind  cc)niponents,  « 
and  r,  resfx-ctively.  Observations  must  be  in  the 
form  of  one  of  these  three  variables  or  in  the 
fonn  of  geopotential  thickness  (the  difference  in 
height  between  two  adjacent  pressure  surfaces). 
For  f'xample.  observations  of  surface  pressure 
are  converted  to  KHK)  mb  heights,  while  satellite 
lemnerature  soundings  are  converted  to  geopo- 
lential  thickness.  Wind  direction  and  .speed 
measurements  are  decomposed  into  u  and  r 
compf)nent,s. 

In  a  data  assimilation  system  like 
N<  ttl.APS,  the  t)-hr  model  forecast  valid  at  the 


analysis  tinte  is  used  as  the  background,  or  first 
guess,  for  the  analysis.  The  analysis,  in  tunt, 
combines  information  from  the  background  with 
the  qucility-controlled  observ'ations  to  provide 
updated  fields  that  are  used  to  initialize  the  fore 
cast  model.  Thus,  the  analysis  is  art 
analysis.  Observation  increments  are  computed 
for  all  data  by  interpolating  the  appropriate 
background  field  to  the  observation  locations 
and  subtracting.  The  irtcremerils  are  then  ana- 
l.vzed  and  added  to  the  background  field  to  pro¬ 
duce  the  latest  three-dimensional  representation 
of  the  atmospheric  height  artd  w  ind  fields. 

The  MVOI  technique  utilizes  data  front 
all  available  platforms  to  produce  artalyzed 
fields  that  possess  minimum  error  in  a  statistical 
sense,  while  taking  into  account  the  cnor  (trop- 
erties  of  the  model  forecast  used  for  the  back¬ 
ground  fields  and  the  error  properties  of  the 
observational  data  used  to  update  the  back¬ 
ground.  The  MVOI  also  constrains  the  results  so 
that  the  analyzed  height  and  wind  increments 
arc  in  geostropic  bahmee  outside  of  the  tropics, 
and  the  vertical  consistency  of  the  model  back¬ 
ground  fields  are  maintained.  Thus,  unlike  earli¬ 
er  objective  analysis  schemes  that  essentially 
averaged  the  nearby  data,  the  MVOI  makes  some 
very  intelligent  decisions  about  how  to  treat  the 
various,  often  conflicting,  sources  of  informa¬ 
tion  it  has  available. 

Four  m^or  factors  influence  the  quali¬ 
ty  of  the  analyzed  product.  One  is  the  quality  of 
the  forecast  used  as  the  analysis  background; 
another  is  the  quality  of  the  observ'ations  them¬ 
selves,  which  requires  the  application  of  sophis¬ 
ticated  quality  coittrol  algorithms.  The  third  fac¬ 
tor  is  the  selection,  during  the  analysis,  of  the 
appropriate  mix  of  mass  and  wind  observations 
from  the  various  sensors.  The  last  factor  is  the 
proper  specification  of  the  statistical  param¬ 
eters  required  by  Uie  optimum  interpolation 
methodology. 

Since  the  skill  of  the  forecast  model  is 
dependent  upon  the  quality  of  the  other  system 
components,  it  is  impossible  to  completely  iso¬ 
late  tmd  define  model  skill.  However,  the  skill  of 
the  data  assimilation  system  ;is  a  whole  is 
reflected  in  Figure  1.  The  -500  mb  anomaly  corre¬ 
lation  is  commonly  used  as  a  measure  of  fore¬ 
cast  skill  by  the  operational  centers.  A  correla¬ 
tion  of  0.6  is  considered  to  be  the  minimum 
value  for  defining  a  useful  forecast.  As  can  be 
seen  from  Figure  1,  the  skill  scores  for  NOGAPS 
show  steady  improv'emeni  over  the  four  year 
p«‘riod  from  1987  to  1990.  It  is  also  interesting  to 
compare  the  relative  skill  of,  for  example,  a 
48-hr  forecast  in  1987  with  a  72-hr  forecast  in 
1990.  Glearly,  there  has  been  a  one-day  improve¬ 
ment,  with  a  three-day  forecast  now  having  the 
same  skill  that  a  two-day  forecast  had  four  years 
ago. 
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FIGURE  1.  NOGAPS  500  mb  height  anomaly  correlation  in  the  Northern  Hemisphere  extratropics  tor  the  years  1987  to  1990. 
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To  ensure  that  erroneous  observations 
are  identified  and  rejecti'd  or  corrected,  data 
quality  control  is  perfonned  in  several  steps. 
First,  except  for  locally  generated  synthetic 
data,  all  observations  are  checked  before  they 
are  passed  to  the  analysis  (Baker,  1991).  These 
tests  include  checks  against  climatological  limits 
and  checks  for  internal  consistency  between 
various  parameters  within  each  report.  For 
example,  if  the  wind  direction  is  reported  as 
variable,  the  wind  speed  must  be  under  7  m  sec  ' 
or  the  wind  observation  is  rejected;  if  both  tem¬ 
perature  and  dew-point  temperature  are  report¬ 
ed,  the  dew-point  temperature  cannot  be  greater 
than  the  air  temperature,  and  over  water,  the  dif¬ 
ference  between  them  is  not  allowed  to  exceed 
30°C.  In  the  case  of  sounding  data,  checks  for 
vertical  consistency  are  also  applied,  including 
lapse-rate  checks,  hydrostatic  checks,  and 
checks  for  excessive  vertical  wind  shears. 
Erroneous  values  are  flagged  for  rejection,  while 
questionable  values  are  flagged  as  suspect. 

The  suspicious  observations  are 
checked  further  within  the  analysis,  after  the 
observations  have  been  converted  to  observed 
increments  by  interpolating  and  subtracting  the 
appropriate  backgrf)und  variable  at  each  obser¬ 
vation  location.  The  tinalysis  gross-error  check 
compares  the  increnu'nts  to  predefined  toler¬ 
ances  that  are  functions  of  the  location  and  type 
of  observ  ation,  the  observation  error  assigtied  to 
the  data,  and  the  iirediction  error  associated 
with  the  background  at  that  location  (Goerss 
and  Phoebus,  1992b).  Increments  exceeding 
these  limits  may  be  njected  outright  or  flagged 
as  suspect.  All  qu»>stionable  increments  are  fur¬ 
ther  checked  against  surrounding  data  during 
the  optimum  interjrolation  analysis.  If  an  incre¬ 
ment  cannot  be  corroborated  by  surrounding 


data,  it  is  excluded  from  the  final  analyzed 
product. 

The  appropriate  mix  of  different  types 
of  observations  from  the  various  sensors  is  guar¬ 
anteed  in  several  ways.  By  selecting  data  and 
performing  the  analysis  within  large  volumes 
that  contain  many  grid  points  (Lorenc,  1981),  we 
avoid  the  problems  associated  with  trying  to 
select  the  few  best  observations  for  each  grid 
point;  we  essentially  allow  all  the  data  over  a 
large  volume  to  influence  the  outcome  at  a  given 
point.  Furthermore,  the  volume  size  is  varied  so 
that  volumes  in  data-sparse  regions  (over  the 
oceans)  are  larger  than  those  in  data-rich  areas 
(over  the  Northern  Hemisphere  continents).  The 
smallest  volumes  typically  cover  an  area  approx¬ 
imately  600  km  on  a  side,  with  the  largest 
around  2,400  km  on  a  side,  and  ail  volumes 
extend  throughout  the  depth  of  the  atmosphere. 
There  are  limits  set  on  the  number  of  observa¬ 
tions  of  one  type  that  will  be  included  in  a  par¬ 
ticular  volume.  Thus,  in  data-rich  areas,  we 
ensure  that  no  one  data  type  wdll  be  selected  at 
the  expense  of  another. 

All  of  the  available  data  within  a  speci¬ 
fied  volume — both  upper-air  and  surface  data — 
are  differenced  from  the  background  field  and 
used  to  compute  the  analyzed  increments  at 
each  grid  point  in  the  volume.  Because  the  anal¬ 
ysis  is  multivariate  and  three-dimensional,  the 
entire  set  of  height,  thickness,  and  wind  incre 
ments  influence  both  the  height  and  wind  fii’ltls. 
The  extent  of  this  influence  is  controlled  uy  the 
various  statistical  parameters  that  are  part  of  the 
optimum  interpolation  formulation 

The  determination  of  the  statistical 
parameters  required  by  the  analysis  is  accom¬ 
plished  using  a  spiecial  database  created  for  the 
purpose  of  monitoring  th''  performance  of  the 
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MVOI.  This  database  contains  a  history  file  of 
observation  minus  background  incren;ents  and 
their  associated  quality  flags,  so  only  increments 
that  passed  all  analysis  quality  che<-ks  are  used 
in  the  calculations.  Since  the  MVOI  is  an  incre¬ 
mental  analysis,  the  required  statistics  must  per¬ 
tain  to  the  itirionenis.  The  statistics  that  must 
be  specified  for  the  MVOI  arc  the  errors  in  the 
background  fields  (prediction  error),  the  predic¬ 
tion  error  <'orrelation  functions  for  all  possible 
combinations  of  data-type  pairs  (e.g.,  hh,  hu  and 
hi;),  and  the  observational  errors  for  each  data 
soui\-e.  The  data  increments  in  the  historical 
database  reflect  both  the  prediction  error  and 
the  observation  error,  but  this  total  error  can  be 
partitioned  into  the  two  separate  error  compo¬ 
nents  following  Hollingsworth  and  Lonnberg 
(1981)).  The  prediction  error  is  normally  deter¬ 
mined  from  the  North  American  radiosonde 
increment  database.  Once  the  prediction  error  is 
known,  the  observation  errors  can  be  estimated 
from  similar  databases  for  each  data  source.  The 
details  of  how  these  quantities  are  derived  are 
discussi'd  in  Goerss  and  Phot'bus  ( 1992b). 

F'or  the  purposes  of  this  paper,  we  note 
that  the  horizontal  and  vertical  prediction  error 
correlation  functions  are  determined  from  actu¬ 
al  data  increments  in  such  a  way  that  the  ana¬ 
lyzed  height  and  wind  increments  will  be 
geostrophically  consistent  outside  of  the  tropics. 
The  prediction  errors  specified  for  height  and 
wind  data  are  also  geostrophically  consistent, 
thus  the  height  errors  vary  somewhat  with  lati¬ 
tude.  Both  height  and  wind  errors  vary  with  vol¬ 
ume  size  to  account  for  the  fact  that  the  back¬ 
ground  field  is  somewhat  less  accurate  in  data 
sparse  areas.  Table  1  shows  sample  values  of  the 
prediction  errors  at  1000  mb,  as  derived  from 
actual  data  increments.  Volume  sizes  range  from 
tiOO  to  2,4(X)  km  along  a  horizontal  side,  depend¬ 
ing  upon  volume  location. 

The  observationtd  error  accounts  for 
both  the  instrument  error  of  Ute  sensor  and  the 
sub-grid  scale  error,  that  is,  how  representative 
the  observation  is  of  the  scales  being  analyzed. 
When  a  new  data  source  becomes  available,  its 
observation  error  must  be  determined  so  that  it 
receives  the  proper  weight  relative  to  the  other 
data  types.  The  weight  a  particular  observation 
receives  at  a  particular  grid  point  is  a  function  of 
its  observ'ation  error,  the  prediction  error  at  the 
observation's  location,  and  the  observation's 
correlation  with  the  grid  point  and  with  the 
other  data  in  the  volume.  For  example,  if  two 
observations  are  co-located,  the  observ'ation 
with  the  larger  observation  error  will  receive  the 
smaller  weight  and  will  have  less  impact  on  the 
outcome  of  the  analyzed  product. 

In  the  past,  the  error  characteristics  of 
the  conventional  marine  surface  observations 
were  treated  as  equal,  mainly  due  to  the  fact  that 


TABLE  1.  Prediction  error  standard  deviations. 


Pressure 

Latitude 

Vol.  Size 

Wind 

'Height 

mb 

m  s  ’ 

m 

1000 

60" 

Small 

20 

105 

1000 

Average 

2.3 

12,0 

1000 

Large 

2.6 

13.6 

1000 

45" 

Small 

2.0 

8  5 

1000 

Average 

23 

98 

1000 

Large 

2.6 

111 

1000 

30" 

Small 

20 

6.0 

1000 

Average 

2.3 

6.9 

1000 

Large 

2.6 

7.8 

the  marine  surface  data  were  stored  and  pro¬ 
cessed  without  regard  to  data  source.  However, 
given  the  importance  of  marine  surface  fore¬ 
casts  to  Navy  operations  and  to  the  success  of 
coupled  air/ocean  forecast  models,  recent 
research  efforts  ha\  e  focused  on  better  specifi¬ 
cation  of  the  error  properties  of  the  marine 
observations  and  of  the  analysis  background 
over  the  oceans. 

Originally,  we  assumed  that  the  height 
and  wind  observ  ation  errors  were  the  same  for 
all  surface  data.  However,  with  tlie  recent  intro¬ 
duction  of  several  new  types  of  observations,  we 
decided  to  revisit  this  problem  and  compute  the 
statistics  separately  for  each  different  marine 
data  source.  Monitoring  the  surface  (observ'ation 
minus  background)  increments  by  data  type  and 
by  region  indicated  that  measurements  made  by 
some  platforms  were  better  than  others.  Thus, 
we  now  identify  the  source  of  each  observation 
and  treat  tlieir  statistics  independently.  The 
changes  in  the  specified  observation  errors  are 
given  in  Table  2,  which  reflects,  for  example, 
that  the  errors  in  the  moored  buoys  reports  are 
smaller  than  those  from  the  ship  reports,  as  are 
the  island  wind  reports.  Elxamination  of  Table  1 
and  Table  2  together  illustrates  how  the  errors  in 
the  model  forecast  compare  to  the  observation 
errors  of  the  data.  The  accuracy  of  global  mod¬ 
els  today  is  such  that  the  observation  errors  in 
the  data  may  actually  exceed  the  prediction 
error  in  the  6-hr  forecast  used  for  the  analysis 
background. 


TABLE  2.  Observation  error  standard  deviations. 

Height  uor  irWind 

m  ms-’ 


Old 

New 

Old 

New 

Island 

60 

7.0 

— 

1.9 

Ship 

60 

70 

2.2 

2.4 

USN  ship 

60 

70 

2.2 

22 

Moored  Buoy 

6.0 

40 

22 

18 

Drifting  Buoy 

6.0 

70 

2  2 

2  2 

C-MAN 

60 

70 

2.2 

22 

PAOB 

24.0 

20.Q 

N.A 

N  A. 

Sic  Bogus 

6.0 

7.0 

2  2 

24 

SSM/I 

N  A 

N  A 

22 

24 
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IRE  2.  Typical  distribution  of  surface  ship,  buoy,  and  C-MAN  stations  within  a  6-hr  data  analysis  window. 
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^INE  DATABASE 

TTiile  tht'  radiosonde  is  the  backbone  of  the 
E  observational  database  over  continenud 
as,  the  number  of  radiosondes  launched  over 
(X'eiuis  is  limited  to  a  few  from  the  more 
ulated  islands  and  an  occasional  ship  of 
ortunity.  These  isolated  upper  air  stations 
supplemented  with  temperature  soundings 
n  the  polar-orbiting  satellites,  cloud-tracked 
d  observations  from  the  geostationary  satel- 
>,  and  aircraft  wind  reports.  However,  at  the 
x-ean  interface,  additional  observations  nrv 
vided  by  military,  research,  and  ctfmmercial 
)s;  moored  and  drifting  buoys;  coastal  marine 
miated  network  (C-MAN)  .stations;  ishmds; 
high-resolution  surface  wind  speed  esli- 
es  from  the  SSM/1.  Because  th(-  MVOl  analy- 
s  a  three-dimensional  analysis,  the  upper-air 
ervations  impact  the  resulting  fields  at  the 
est  analysis  level,  but  this  influence  is  limited 
irface  data  are  available  in  the  same  area, 
s,  we  will  only  describe  the  surface  data  in 
paper.  Additional  infonnation  about  the 
re  database  is  givmi  by  Coerss  and  Phoebus 
»2a). 

A  marine  surface  observation  may  be 
isured  or  considered  valid  at  heights  that 
I  from  a  few  to  several  tens  of  meters  above 
actual  air/ocean  interface,  depending  on  the 
>  and  location  of  the  instrument  making  the 
Lsuroment.  Since  the  vertical  resolution  of 
global  model  does  not  support  these  distinc 
s,  all  marine  surface  observations  are  con¬ 
ed  to  observations  of  lOdO  mb  height  and 
T  mb  u  and  c  wind  components.  We  do, 

'ever,  make  u.se  of  boundary  layer  informa- 
when  we  do  this  assignment.  The  model 


produces  a  10  m  wind  field  that  is  stability 
dependent  and  consistent  with  the  model’s 
boundary  layer  physics,  and  we  have  found  that 
these  winds  agree  more  closely  with  tlie  marine 
data  than  the  1000  mb  winds  do.  Thus,  to  con- 
\  ert  the  surface  wind  measurements  to  1000  mb 
increments,  we  actually  compute  boundary  layer 
wind  increments  by  interpolating  and  subtract¬ 
ing  the  model’s  10  m  winds  from  the  observed 
surface  wind  components.  These  10  m  u  and  i> 
wind  increments  are  then  assigned  to  the  1000 
tub  level.  By  assuming  that  the  wind  increments 
do  not  \  ary  vertically  over  the  small  distance 
between  10  m  and  1000  mb,  we  are  preserving 
the  vertical  structure  defined  by  the  model’s 
boundary  layer. 

Conventional  Data 

The  conventional  ntarine  surface  data 
include  reports  from  ships,  moored  and  drifting 
buoys,  and  C-MAN  stations.  A  typical  1200  uni¬ 
versal  time  clock  (IITC)  analysis  product  might 
include  400  or  5(K)  ships,  slightly  less  than  100 
moon'd  and  drifting  buoys  (each),  and  50  or  60 
C-MAN  reports.  Figure  2  illustrates  the  global 
distribution  of  these  observations. 

Approximately  90  percent  of  the  reports  are  in 
the  Northern  Hemisphere.  The  variables  normal¬ 
ly  reported  by  these  marine  platforms  include 
sea-level  pressure,  temperature,  moisture,  wind 
speed,  and  wind  direction.  Drifting  buoys  are  the 
exception  since,  for  the  most  part,  the  only 
atmospheric  vairiable  they  measure  is  sea-level 
pressure. 

Many  of  the  marine  platforms  report 
frequently.  As  the  data  arc  processed,  duplicate 
or  redundant  reports  are  identified  by 
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FIGURE  3.  Distribution  of  wind  observations  by  wind  speed 
category.  Percentage  of  SSM/I  observations  in  each  category 
compared  to  percentage  of  ship  and  buoy  observations  in 
the  same  categories.  _ 


WIND  SPEED  CATEGORY  (m/s) 

comiraring  latitudes,  longitudes,  and  hlock/sta- 
I'on  names.  Reports  at  the  same  location  or  with 
the  same  identifier  are  removed.  Only  one  obser¬ 
vation  per  platfonn  is  used  in  the  analysis,  even 
if  that  platfonn  reports  every  hour.  The  observa¬ 
tion  saved  is  the  one  closest  to  the  analysis  time. 

Island  Data 

Like  other  surface  observations  over 
land,  islands  typically  report  wind  speed  and 
direction  and  either  sea-level  pressure,  station 
pressure,  or  geopotential  height  at  a  standard 
pressure  surface.  If  pressure  is  reported,  it  is 
converted  to  a  height  (Goerss  and  Phoebus, 
1992b).  Zonal  and  meridional  wind  components 
are  computed  from  the  reported  wind  speed  and 
direction.  Until  recently,  surface  w'ind  observa¬ 
tions  over  land  were  not  used  in  the  analysis. 
Such  observations  are  often  not  representative 
of  the  wind  in  the  free  atmosphere  because  of 
the  effects  of  terrain  or  shallow  radiation  inver¬ 
sions.  However,  many  of  the  island  instrument 
stations  have  been  relocated  to  reduce  these 
(iroblems.  Personnel  from  the  Joint  Typhoon 
Wanting  Center  (,fTWC)  in  Guam,  who  routinely 
monitor  and  analyze  the  Pacific  island  wind 
reports,  indicated  that  .JTWC  generally  accepts 
the  island  wind  directions  as  accurate,  but 
makes  corrections  to  the  reported  speed.  So, 
these  observations  were  re-evaluated  for  use  in 
N(  IGAJ’S.  The  observation  minus  background 
differences  indicated  that  the  Pacific  island  wind 
observations  were  as  accurate  as  the  ship  wind 
reports.  Thus,  these  particular  island  wind 
observations  are  now  included  in  NOGAPS,  but 
only  if  the  reported  wind  speed  is  greater  than  4 
m/s.  This  caveat  rejects  the  light  and  variable  or 
calm  w  ind  reports  that  are  so  common  during 


the  night  but  includes  reports  that  are  more  r('[)- 
resentative  of  a  well-mixed  boundary  layer. 

The  Pacific  islands  cover  an  area  that 
extends  from  3()°S  to  3(PN  and  includes  slightly 
less  than  UK)  island  .stations.  Islands  outside  this 
region  are  still  excluded,  since  they  are  more 
isolated  and  it  is  extremely  difficult  to  delennine 
the  quality  of  their  ob.servations.  Surface  w  ind 
reports  from  continental  stations  are  also 
ignored. 

SSM/I  Wind  Data 

FNOU  uses  the  algorithm  of 
Goodlx'flet  et  al.  ( 1990)  to  estitnate  the  ocean 
surface  (19..')  m)  wind  speeds  from  the  SS,\1/1 
measured  brightness  temperatures.  These 
speeds  range  up  to  30  m  s  '  and  meet  the  siteci- 
fied  accuracy  of  ±2  m  s  '  in  rain-free  conditions. 
Since  heavy  concentrations  of  water  \  apor  ;uid 
rain  effectively  mask  microwave  emission  from 
the  oceiui  surface  at  the  SSM/1  frequencies,  wind 
speed  retrievals  under  these  conditions  are  more 
likely  in  error.  FNOC  uses  four  rain-flags  to  indi¬ 
cate  accuracy,  ranging  from  0  (no  rain,  error  <  2 
m  S“')  to  3  (heavy  rain,  error  >  10  m  s  '). 
Currently,  the  NOGAPS  data  pre-processor 
screens  the  SSM/I  data  so  that  only  wind  obser¬ 
vations  with  a  rain-flag  of  0  or  1  are  used  in  the 
M\'OI. 

The  SSM/1  data  are  valued  not  only  for 
their  accuracy  but  also  for  their  distribution. 
Figure  3  shows  the  distribution  of  the  SSM/1 
wind  speeds  compared  to  those  from  the  ships 
and  buoys.  The  m^ority  of  the  SSM/1  wind 
speed  data  (>  50  percent)  fall  in  the  5  to  10  m  s-‘ 
category,  while  the  0  to  5  m  s-'  category  has  the 
largest  percentage  of  ship  and  buoy  observa¬ 
tions.  What  is  espet  ially  noticeable  is  that,  over¬ 
all,  the  higher  w  ind  speed  categories  contain  a 
larger  percentage  of  SSM/I  observations  than 
ship  observations.  Since  ships  are  deliberately 
routed  away  from  high  wind  speed  areas,  this 
difference  most  likely  reflects  the  indiscriminate 
global  coverage  of  the  SSM/1  data. 

Reasonable  global  coverage  is  provided 
by  the  SSM/1  within  a  twenty-four  hour  period, 
thus  providing  an  important  supplement  to  the 
other  relatively  sparse  observations  over  the 
ocean.  On  the  other  hand,  a  tyT)ical  six-hour  cov¬ 
erage  (the  data  assimilation  cycle)  produces  in 
exce.ss  of  200,000  observations  at  approximately 
2.5  km  resolution,  a  resolution  much  higher  than 
the  1 .5°  resolution  of  NOGAPS.  Tlius,  it  is  rea¬ 
sonable  to  filter  the  information  provided  by  the 
SSM/I  observations.  Since  the  SSM/I  observa¬ 
tions  arc  extremely  consi.stent,  only  eveiy  sixth 
SSM/I  observ’ation  is  pa.ssed  to  the  analysis. 
Furthermore,  any  wind  speed  reported  outside 
the  range  of  the  SSM/1  instrument  is  ignored. 

Next,  the  wind  speed  data  are  subjecl- 
('d  to  quality  control  procedures  that  compare 
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ABLE  3.  SSl^/l  quality  control. 

Preliminary  Rejection  Percentage 

Analysis  Criteria  Rejected 

cWms’  SSM/I>7ms-'  -  5;o% 

1-10  m  s  *  I  Oiff  1>  7.5  ms-'  -  1 .5% 

!lOms-'  I0i«l>10ms-'  -0.1% 


ihe  obson  ed  wind  speed  to  a  preliminary  analy¬ 
sis  of  wind  speed  valid  at  10  in.  The  10  m  wind 
analysis  product  is  obtained  from  a  preliminary 
analysis  that  is  executed  from  the  .same  software 
and  the  same  background  fields  as  the  NOGAPS 
analysis,  but  it  is  executed  at  only  two  hours 
past  the  synoptic  time,  so  there  are  fewer  obser¬ 
vations  available  than  there  are  during  the  nor¬ 
mal  assimilation  cycle.  The  preliminary  analysis 
u.ses  almost  all  of  the  available  data,  except  for 
aircraft  wind  reports,  satellite  temperature 
soundings,  and  the  SSM/I  wind  speed  data,  and 
it  is  perfonned  only  at  the  lowest  four  analysis 
levels  ( 1000  to  700  mb).  By  the  same  reasoning 
that  is  used  to  assign  surface  data  increments  to 
1000  mb,  we  create  the  10  m  wind  analysis  prod¬ 
uct  by  adding  the  analyzed  1000  mb  increments 
to  the  10  m  background  fields. 

Wdiether  or  not  an  observation  is 
rejected  during  the  comparison  to  the  prelimi¬ 
nary  10  rn  wind  product  depends  upon  both  the 
observed  wind  speed  and  the  analyzed  wind 
speed  valid  at  the  observation  location.  The 
rejection  criteria  are  given  in  Table  0,  where 
I  Uiff  I  represents  the  absolute  value  of  the 
observed  minus  analyzed  wind  speed,  and  the 
percentage  rejected  is  in  reference  to  the  total 
number  of  SSNl/l  observations  processed.  The 
strictest  criteria  are  applied  and  the  most  data 


are  rttjected  when  the  analyzeti  w  ind  speed  is 
very  light  coupled  with  an  SSM/1  wind  speed 
that  is  substantially  stronger.  Since  the  prelimi¬ 
nary  analysis  product  will  also  be  used  to  assign 
wind  directions  to  the  SSM/I  observations,  this 
decision  reflects  the  fact  that  it  would  be  inap¬ 
propriate  to  assign  a  wind  direction  from  a  light 
and  variable  condition  to  a  stronger  wind  speed. 
Doing  so  could  create  unrealistic  low-le\’e!  wind 
shears  and  convergence/divergence  zones.  The 
other  categories  of  comparison  exhibit  more  tol¬ 
erance  for  differences  between  the  preliminary 
analysis  fields  and  the  SSM/1  data  and,  as  a 
result,  reject  fewer  observations. 

(omputation  of  Superobs 

Since  the  small  scales  sampled  by  the 
■SS.M/l  are  of  little  importance  to  a  global  data 
assimilation  system,  it  is  reasonable  to  further 
filter  the  information  provided  by  the  SSM/I 
observations.  To  reduce  the  number  of  data 
points  without  losing  infonnaiion  on  scales  rele¬ 
vant  to  NOGAPS,  the  SSM/1  data  are  combined 
\nio  superobs,  which  are  essentially  local  aver¬ 
ages  of  groups  of  observations.  A  typical  distri¬ 
bution  for  the  superobs  available  for  a  single 
analysis  product  is  shown  in  Figure  4,  which 
w  as  plotted  using  a  six -hour  data  window  cen¬ 
tered  at  1200 IITC. 

Before  the  superobs  are  computed,  the 
10  m  w  ind  speeds  from  the  preliminai>  analysis 
fields  are  interpolated  to  the  location  of  each 
individual  observation,  and  the  interpolated  val¬ 
ues  are  subtracted  from  the  observed  wind 
speeds  to  form  SSM/I  wind  speed  increments. 
Nexd,  the  globe  is  divided  into  200  km  boxes, 

;ind  the  SSM/1  wind  speed  iucremctils  within 


FIGURE  4.  Typical  distribution  of  SSM/I  superobs  within  a  6-hr  data  analysis  window. 
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each  box  are  averaged.  The  increments  rather 
than  the  original  values  are  used  becaust»,  in 
general,  tlie  increment  field  is  more  horizontally 
homogeneous  than  the  raw  data.  Once  the 
superob  increments  are  computed,  they  are 
assumed  to  be  valid  at  the  centers  of  the  boxes. 
Finally,  the  superob  wind  speeds  are  detennined 
by  adding  each  superob  increment  to  tlie  10  ni 
analyzed  wind  speed  interpolated  to  the  superob 
location. 

Statistical  comparisons  of  obst'rved 
minus  analyzed  increments  for  the  original 
SSM/1  wind  speed  data  and  the  SSM/1  wind 
speed  superobs  (Table  4)  illustrate  that  the  dif¬ 
ferences  in  the  root-mean-square  (rms)  errors 
are  small,  which  indicates  that  the  superob  pro¬ 
cess  is  not  introducing  significant  errors  in  the 
data.  If  anything,  the  averaging  proce.ss  tends  to 
reduce  the  random  noise  that  is  present  in  indi¬ 
vidual  observations.  Sintilar  statistics  were  pro¬ 
duced  for  superobs  classified  by  the  number  of 
indiridual  SSM/I  observations  composing  the 
superobs.  The  errors  were  larger  for  those  cases 
where  only  a  few  increments  were  averaged — 
most  likely  near  coastlines  or  along  the  edges  of 
swaths,  where  the  SSM/I  data  may  not  be  as  reli¬ 
able.  Furthermore,  the  error  introduced  by 
assigning  the  superob  to  a  location  at  tfte  center 
of  the  box  will  be  greater  when  there  are  fewer 
daua  points  in  the  box.  Thus,  superobs  are 
computed  only  for  boxes  with  at  least  five 
observations. 

Assignment  of  Wind  Direction 

The  analysis  variables  are  geopotential 
height  and  the  u  and  v  wind  components — vari¬ 
ables  that  have  well-known  dynamical  relation¬ 
ships.  Thus,  the  MVOI  analysis  does  not  utilize 
scalar  wind  speed  information  directly.  So.  once 
the  SSM/1  superob  increments  have  been  com¬ 
puted,  each  superob  wind  speed  must  be 
assigned  a  wind  direction.  There  have  been  sev- 


TABLE  4.  SSM/I  Wind  speed  statistics.  May  24-28, 1991. 

Original  SBM/l  Wind  Speeds 

Ob  —  Prelim^  Anal,  (n^r'l 


SSM/I 

(ms') 

Mean 

St  Oev 

rms 

0-  5 

-12 

1.8 

24 

5-10 

04 

1.7 

1.7 

10-15 

1.9 

2.1 

28 

15-20 

36 

2  1 

4.2 

20-25 

6.4 

22 

6.9 

Superob 

SSM/1  Superob  Wind  Speeds 

Ob  —  Prelim.  Anal,  (m  s-’j 

(m  s') 

Mean 

St  Oev 

rms 

0-  5 

-0.4 

1.5 

1.5 

5-10 

06 

15 

16 

10-15 

1.9 

20 

2  7 

15-20 

35 

1.9 

4  0 

20-25 

6.0 

2.3 

6.3 

eral  techniques  profxxsed  for  indr-jKmdently 
detennining  appropriate  directions  for  the  SSM/1 
wind  speeds.  Atlas  et  al.  (1991 )  combine  conven¬ 
tional  wind  data  and  SSM/I  winds  with  back¬ 
ground  tmalyzed  surface  wind  fields  from  the 
Fluropean  (tentre  for  Medium  Range  Weather 
F'orecasLs,  using  variational  teeiiniques  to  con¬ 
strain  the  degri'e  to  which  the  SSM/1  data  can 
modify  the  background.  Yu  ( 1SI87)  bases  his 
technique  on  Kknian  boumlary  layer  dynamics, 
using  the  National  Meteorological  Center's  sur¬ 
face  pressure  analysis  product  and  the  SSM/I 
wind  speed  to  compute  the  surface  drag  coeffi¬ 
cient,  which  then  allows  the  computation  of  the 
vector  wind  at  some  height  z*  within  the  bound¬ 
ary  layer  (tyjtically  10  m). 

In  our  case,  we  use  the  pow'cr  of  a 
global  data  assimilation  system  to  produce  the 
10  m  preliminary  fields  using  the  most  recent 
data,  independent  of  the  SSM/1.  This  tmalyzed 
wind  field  is  inleri)olated  to  the  location  of  each 
superob  wind  speed,  and  the  resulting  wind 
directions  are  assigned  to  the  superob  locations, 
FYom  this  point  on,  the  SSM/I  superob  compo¬ 
nents  are  treated  like  any  other  surface  wind 
observations. 

S3Tithetic  Data 

Final  y,  the  global  database  is  supple¬ 
mented  by  the  creation  of  synthetic  observations 
for  specific  applications.  The  most  commonly 
used  syntltetic  observations,  which  are  available 
to  all  the  operational  forecast  centers,  are  the 
subjectively  derived  estimates  of  sea-level  pres¬ 
sure  (PAOBS)  that  are  produced  by  the 
Australian  Bureau  of  Meteorology  from  satellite 
imagery.  These  observations  are  evenly  dis¬ 
tributed  over  both  land  and  sea,  but  are  confined 
to  the  Southern  Hemisphere.  NOGAPS  uses  this 
data  only  south  of  20°S.  Approximately  150 
PAOBS  are  available  at  FNOC  for  the  0000  UTC 
and  1200  UTC  analysis  products. 

Locally,  FNOC  personnel  generate  syn¬ 
thetic,  or  bogus,  surface  observ'ations  in  the 
vicinity  of  oceanic  extra-tropical  cyclones. 
Typically,  such  observations  are  created  in  rela¬ 
tively  data-sparse  areas  like  the  North  Pacific,  in 
situations  where  satellite  imagery  indicates  that 
the  cyclone  is  deeper  than  its  depiction  by 
NOGAPS.  This  technique  is  particularly  useful 
for  capturing  explosively  deepening  lows  in  the 
Gulf  of  Alaska,  since  in  the  absence  of  other 
data,  the  synthetic  observations  can  provide  a 
surface  baseline  for  the  layer  thicknesses 
derived  from  the  satellite  soundings. 

Only  the  central  location  and  sea-level 
pressure  of  the  low  is  subjectively  specified.  The 
difference  between  this  synthetic  observation 
and  the  analyzed  sea-level  pressure  at  the  same 
location  is  then  used  to  automatically  generate 
eight  additional  synthetic  observations,  centrally 


located  325  Ion  from  the  storm’s  center  on  radii 
separated  by  45°  (Goerss,  1989).  These  outer 
observations  include  both  sea-level  pressures 
and  winds.  Tlie  pressures  are  determined  by 
assuming  that  the  pressure  difference  changes 
unifonnly  from  the  central  difference  to  zero  as 
distance  from  the  center  increases  from  0  to  650 
km  along  a  given  radius.  Thus,  half  the  central 
pressure  difference  is  added  to  the  analyzed  sea- 
level  pressure  at  each  of  the  eight  observation 
locations.  The  sjnthetic  wind  observations  are 
determined  by  adding  the  gcostrophic  compo¬ 
nents  implied  by  the  pn?ssure  difference  distri¬ 
bution  to  the  antdyzed  wind  components  at  the 
observation  locations. 

Once  the  synthetic  surface  observa¬ 
tions  have  been  generated,  they  are  temporarily 
entered  into  the  NOOAl’S  operational  database. 
A  preliminan,'  imalysis  is  then  perfonued  with 
the  MVOI  for  a  limited  area  surrounding  the 
cyclone  of  interest.  If  the  P'NOC  personnel  are 
satisfied  with  the  resulting  changes  in  the  low, 
the  synthetic  observations  are  entered  into  the 
dabilutse  permanently,  and  thcj  are  accessed  for 
the  iK'xt  analysis/forecast  update  cycle. 

After  changes  to  the  physical  parame- 
terizations  and  resolution  of  the  forecast  model 
in  1989,  NOGAPS  began  actively  producitig  and 
developing  tropical  circulations.  To  ensure  that 
the  circulations  were  maintained  in  the  proper 
location  and  at  the  corrc'ct  strength  throughout 
the  data  assimilation  cycle,  NRL  developed  a 
technique  for  assimilating  synthetic  v'hid 
soinidiiigs  in  the  vicittity  of  these  tropical 
stonns  (Goerss  et  al.,  li)91).  In  this  case,  the 
generation  of  the  synthetic  observations  is  cont- 
pletoly  automated,  using  the-waniing  messages 
i.ssued  by  .JTW((  and  by  the  National  Hurricane 


FIGURE  5.  Location  of  synthetic  tropical  soundings  relative 
to  storm  center,  overlayed  on  an  850  mb  height  analysis  of 
typhoon  Pat. 
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Center  in  Miami  for  tropical  cyclones  whose 
maximum  wind  speeds  are  greater  than  13  m/s 
(depression  stage).  The  synthetic  soundings  arc 
centered  on  the  stonn’s  location  and  reflect  the 
maximum  wind  spo<‘d  and  radius  of  maxiimmi 
wind  information  contained  in  the  warning  mes¬ 
sage. 

Each  synthetic  tropical  sounding  con¬ 
tains  wind  infonnation  from  1000  mb  to  4(M)  mb, 
along  with  a  1000  mb  height  observation.  If  the 
circulation  has  reached  tropical  storm  strength 
or  greater,  soundings  are  generated  at  thirteen 
points;  one  at  the  center;  four  located  220  km 
north,  south,  east,  and  west  of  the  center;  four 
located  440  km  northeast,  southeast,  southwest, 
and  northwest  of  the  center;  and  four  locat«'d 
660  km  north,  south,  cast,  and  west  of  the  center 
(Figure  5).  For  tropical  depressions,  soundings 
lire  generated  only  at  the  five  most  interior 
points.  All  positions  are  updateti  and  new  olisei - 
\'alions  are  generated  every  6  hrs. 

Each  sounding  is  composed  of  two 
parts.  The  first  is  the  contribution  of  the  large- 
scale  environmental  flow,  determined  by  spec¬ 
trally  truncating  the  NOGAPS  background  field 
.so  that  only  the  first  twenty  wave  numbers  are 
retained.  The  other  is  the  vortex  comj)onent, 
derived  from  a  symmetric  Rankine  vortex  in  gra¬ 
dient  balance  (Anueisson  and  Hollingsworth. 
19nS).  The  structure  of  the  Rankine  vortex  is 
controlled  by  the  maximum  wind  speed,  the 
radius  of  maximum  w'md,  and  an  ex^ponential 
factor  the  governs  the  flow  beyond  this  radius. 
The  contribution  of  the  symmetric  vortex  com¬ 
ponent  is  reduced  in  the  vertical  to  simulate  the 
warm-core  structure  of  the  tropical  cyclone. 
Furthermore,  the  Rankine  vortex  wind  speeds 
do  not  drop  off  quickly  enough  with  distance  in 
the  horizontal.  Thus,  the  vortex  wind  speeds  in 
the  outer  two  circles  are  reduced  by  40  percent 
at  (KiO  km  and  by  20  percent  at  440  km  to  com¬ 
pensate.  The  resulting  set  of  synthetic  observa¬ 
tions  produces  a  tropical  cyclone  that  is  consis- 
tt  (It  with  how  the  global  model  depicts  such  fea¬ 
tures  rather  than  how  tropical  storms  actually 
appear  in  nature. 

DATA  IMPACT 

Impact  of  SSM/I  Data 

Tlte  performance  of  the  data  assimila¬ 
tion  .systeui  is  monitored  by  continually  produc¬ 
ing  statistics  of  observation  minus  background 
and  observation  minus  analysis  differences.  Just 
after  the  introduction  of  the  SSM/I  data  into 
NOGAPS,  global  statistics  were  computeft  from 
the  observation  minus  (KXK)  mb)  background 
differences  (Table  5).  First,  we  notice  that  a  sig¬ 
nificant  bias  is  present  at  the  higher  wind 
speeds,  which  is  also  reflected  in  the  rms  error. 
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TABLE  5.  SSIWI  wind  speed  statistics.  September  13-14,  iggo. 


SSM/1 
(m  s’) 

Mean 

St  Dev 

rms 

0-  5 

-1.7 

22 

28 

5-10 

-0.3 

23 

23 

10-15 

0.8 

33 

3.4 

15-20 

2.7 

3.1 

4.1 

20-25 

66 

24 

7.0 

TABLE  6.  SSM/I  wind  speed  statistics.  September  23-24, 1990. 


SSM/I 
(111  s  ') 

Mean 

St  Dev 

rms 

0-  5 

-15 

20 

25 

5-10 

-0.2 

2  1 

21 

10-15 

08 

27 

28 

15-20 

1  8 

26 

32 

20-25 

4.3 

26 

50 

clilToreiift's  have  decreased  over  lime.  Thus,  we 
conclude  that  the  forecast  system  is  aceeitling 
the  data.  As  a  means  of  indep(*ndent  veririealioii, 
we  also  computed  the  vector  nus  ('rrors  ol 
observation  minus  background  ineremt'nts  for 
the  ship  and  buoy  data  before  and  iifler  the 
inclusion  of  the  SSM/1  wimi  speeds  From  'fabU' 
7,  it  is  apparent  that  the  ( 1000  mb)  background 
also  agrees  more  closely  with  these  independent 
data  .sources  after  the  SSM/I  data  were  intro 
duced.  The  largest  enange  is  observed  in  the 
Southeni  Hemisphere,  as  would  be  expected, 
since  any  new  data  set  will  haw  a  more  signifi. 
cant  impact  in  data  sparse  art'as. 

To  make  the  comparison  to  ifidepen 
dent  data  sources  more  consistt'nl,  we  cone 
pared  the  SSM/1  superob  wind  speed  smtislies  to 


TABLE  7.  Ship/buoy  wind  statistics. 

0*'  —  Background  Differences 
rms  Before  SSM/I  Alter  SSM/1 

0  Hem  5.0  ms-'  4.9  ms-' 

SHem  5,5  ms-'  5.1ms-' 


TABLE  8.  Marine  wind  speed  statistics  May  24-28, 1991. 

SsM/l  $uperob  Wind  Speeds 

Ob  —  Background  (m  s-') 


Ship/Buoy  Wind  Speeds 

Ob  —  Background  (m  s-') 


FIGURE  6.  Difference  fields  of  (a)  100P  mb  wind  vectors 
(maximum  vector  corresponds  to  5  m/s)  and  (0)  sea-level 
pressure,  contoured  e  y  mb.  Differences  are  for  an 
analysis  with  SSM/I  data  minus  an  analysis  without  SSM/I 
wind  data. 
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How  ever,  a  similar  bias  is  present  in  the  statis¬ 
tics  for  other  ty{)es  of  surface  w  ind  data  incre¬ 
ments,  indicating  that  the  model  tends  to  under¬ 
estimate  the  .sf longer  wind  speeds.  This  bias  is 
most  likely  due  to  the  fact  that  mary  of  these 
high  wind  speed  events  cannot  be  well  repre¬ 
sented  by  the  resolution  of  the  current  model. 
However,  the  source  of  this  bias  is  under 
investigation. 

Ten  days  after  NOGAPS  incorporated 
the  SSM/I  data,  the  sante  statistical  computa¬ 
tions  were  performed  once  again  (Table  6). 
Comparison  to  Table  r>  shows  that,  for  all  wind 
speeds,  the  observation  minus  background 
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{ 

_  the  same  statistics  for  ship  and  buoy  wind 
I  speeds  for  the  1200  UTC  analyses  from  May  24 

«  to  28,  1991.  The  mean,  standard  deviation,  and 

mis  error  for  the  observ'ation  minus  (10  m) 
background  differences  are  shown  in  Table  8. 
The  SSM/I  data  do  not  appear  to  be  biased  rela¬ 
tive  to  the  other  marine  surface  wind  observa¬ 
tions,  and  their  difference  statistics  are  compa¬ 
rable.  WTien  compared  to  Table  4,  it  is  also  clear 
that  the  SSM/1  wind  speeds  agree  more  closely 
with  the  preliminary  analysis  than  with  the  anal¬ 
ysis  background  field,  further  showing  consis¬ 
tency  between  the  SSM/I  data  and  the  other 
marine  wind  observations.  While  the  statistics 
presented  here  are  from  a  limited  time  period, 
there  are  only  slight  day-to-day  variations  in 
these  global  values. 


FIGURE  7.  Difference  fields  of  (a)  1000  mb  wind  vectors  (maximum  vector  corresponds  to  10  m/s)  and  (b)  1000  mb  heights, 
contoured  every  10  m.  Differences  are  for  a  48-br  forecast  made  iv/f/r  minus  a  48-hr  forecast  made  without  SSM\  wind  data 
In  the  Without  SSM/I  data  case,  the  observations  were  withheld  for  48  hrs  of  the  data  assimilation  cycle  prior  to  making  the 
comparison  forecast. 


Finally,  Figure  6  shows  the  impact  of 
the  SSM/I  wind  speeds  on  the  analysis  of  a  low 
pressure  system  in  the  Southern  Hemisphere. 
Inclusion  of  the  SSM/I  data  increases  the 
strength  of  the  flow  around  the  low  (Figure  6a), 
and  through  the  multivariate  coupling  of  the 
analysis,  consistently  decreases  the  central  pres¬ 
sure  of  the  low  by  2  mb  (Figure  6b).  While  this 
may  seem  like  a  small  difference,  such  differ¬ 
ences  in  the  analyzed  fields  can  often  make  large 
differences  in  the  model  forecasts.  As  an  exam¬ 
ple,  SSM/I  data  were  withheld  from  the  data 
assimilation  cycle  for  48  hrs;  then  a  48-hr  fore¬ 
cast  was  made  from  these  initial  fields.  The  fore¬ 
cast  without  the  SSM/1  data  was  differenced 
from  the  NOGAPS  48-hr  forecast  valid  at  the 
same  time  from  the  operational  nm.  As  can  be 
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FIGURE  8.  NOGAPS  depiction  ot  supertyphoon  Page  prior  to 
the  assimilation  of  any  s,nthetic  observations.  The  small 
hurricane  symbol  marks  the  actual  storm  position  (a) 
Analysis  valid  at  0000  UTC,  November  19, 1990;  (b)  72-hr 
forecast  valid  at  0000  UTC,  November  22, 1990;  (c)  120-hr 
forecast  valid  at  0000  UTC,  November  24, 1990. 


seen  from  Figure  7,  the  differences  in  I  lie  wind 
field  are  as  large  as  20  ni  s  ‘  and  the  1000  mb 
height  field  exhibits  differences  as  large  as  H2  in. 

This  height  difference  would  translate  to  about 
an  1 1  mb  difference  in  surface  pressure,  just 
from  excluding  SSM/1  data  for  48  hrs.  Although 
the  impact  of  the  SSM/I  data  appears  to  hr-  the 
greatest  near  the  poles  and  in  the  northem  mid¬ 
latitudes,  sizeable  differences  are  also  seen  in 
the  tropical  winds  in  the  western  Pacific  and 
Indian  oceans. 

Impact  of  Synthetic  Tropical  Soundings 

In  1989,  the  resolution  of  the  forecast 
model  was  increased  and  improvements  were 
made  in  the  parameterization  of  the  diabatic 
processes  (Hogan  and  Rosmond,  1991).  As  a 
result  of  these  changes,  NOGAI’S  began  actively  i 

developing  and  forecasting  tropical  stonns.  The 
ability  of  the  forecast  model  to  develop,  intensi 
fy,  and  forecast  the  movement  of  tropical  circu¬ 
lations  accentuated  the  need  for  properly  [losi-  | 

tioning  the  storms  in  the  NOGAPS  initial  fields.  ' 

It  was  also  a  necessary  ingredient  for  successful 
tropical  data  assimilation.  If  the  model  could 
maintaun  these  initial  circulations  throughout  j 

the  forecast  period,  then  synthetic  observations  | 

could  bo  used  to  insert  the  developing  stonns  j 

into  the  analyzed  fields. 

To  illustrate  the  model’s  capabilities 
prior  to  the  use  of  synthetic  data,  Figure  8  dis¬ 
plays  the  NOGAPS  sta-level  pressure  analysis 
product  and  the  72-hr  .md  120-hr  forecast  prod¬ 
ucts  made  from,  this  -analysis.  At  0000  UTC, 

November  19, 1990,  a  1006  mb  low  is  analyzed 
about  500  km  from  the  JTWC  position  of  a  tropi¬ 
cal  depression  that  eventually  developed  into 
supertyphoon  Page  (winds  exceeding  66  m/s). 

The  72-hr  NOGAPS  forecast  valid  at  0000  UTC 
on  November  22  has  developed  a  999  mb  low. 

By  this  time,  in  reality  Page  had  reached  tropica) 
storm  strength  with  maximum  winds  of  18  m/s, 
and  the  error  in  the  forecast  position  is  approxi¬ 
mately  750  km.  After  120  hrs,  the  storm  has 
deepened  to  a  993  mb  low  (Figure  8c).  At  0000 
UTC,  November  24,  Page  was  actually  located 
about  550  km  from  the  forecast  position  tmd  had 
maximum  wind  speeds  of  31  m/s.  Even  though 
the  error  in  the  forecast  position  is  fairly  large, 
the  model’s  development  of  the  storm  is  notable. 

Consider  the  model’s  depiction  of  the 
same  storm  a  few  days  later  (Figure  9).  In  this 
case,  synthetic  observations  have  been  gener¬ 
ated  for  Page  since  it  reached  tropical  storm 
strength  at  0000  UTC,  November  22.  These 
observations  were  assimilated  into  NOGAPS 
beginning  at  1200  UTC,  November  22.  The  0000 
UTC  analysis  product  on  November  24  shows  a 
1002  mb  low  located  approximately  1(X)  km  from 
the  JTWC  position  of  Page,  which  at  this  time  is 
approaching  typhoon  strength  with  maximum 
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FIGURE  9.  NOGAPS  depiction  of  supertyphoon  Page  after 
the  assimilation  of  synthetic  observations.  The  small 
hurricane  symbol  marks  the  actual  storm  position,  (a) 
Analysis  valid  at  0000  UTC,  November  24, 1990;  (b)  72-hr 
forecast  valid  at  0000  UTC,  November  27, 1990;  (c)  120-hr 
forecast  valid  at  0000  UTC,  November  29, 1990. 
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wind  speeds  of  31  iii/s.  The  72-hr  forecast  devel¬ 
ops  a  985  mb  low  located  about  300  kin  from  the 
actual  storm  position.  By  this  time,  0000  UTC, 
November  27,  Page  was  at  its  peak  intensity 
with  maximum  winds  of  72  m/s.  The  120-hr  fore¬ 
cast  shows  a  973  mb  low  whose  center  is  locat¬ 
ed  about  400  km  from  the  verifying  position  of 
Page  at  0000  UTC,  November  29.  The  stonn  was 
still  a  damaging  typhoon  as  it  made  landfall  in 
Japan  on  November  30. 

Comparing  Figures  8  and  9,  the  analy¬ 
sis/forecast  products  with  no  synthetic  data  had 
position  errors  of  500  km,  750  km,  and  550  km, 
while  the  later  senes  had  track  errors  of  100  km, 
300  km,  and  400  km.  While  the  initial  fields  of 
surface  pressure  differed  by  only  4  mb,  the 
series  of  forecast  products  generated  from  the 
analysis  using  synthetic  data  produced  a  storm 
whose  central  pressure  was  20  mb  lower  at  120 
hrs.  Thus,  the  assimilation  of  the  synthetic 
obser\ations  produced  a  more  accurate  repre¬ 
sentation  of  the  stonn  in  the  analyzed  field, 
resulting  in  more  accurate  forecasts  of  the 
.storm’s  track  and  its  strength.  Even  though  the 
circulations  generated  by  NOUAPS  are  larger 
than  those  observed  in  nature,  clearly  the  fore¬ 
cast  model  develops  and  intensifies  tropical  cir¬ 
culations.  This  case  illustrates  a  very  dramatic 
deepening  that  is  not  typical  of  most  tropical  cir¬ 
culations.  However,  in  the  few'  cases  where  such 
a  strong  deepening  has  been  observ'ed  in  the 
NOGAPS  forecasts  prior  to  the  assimilation  of 
synthetic  data,  almost  all  were  associated  with 
storms  that  approached  or  reached  super¬ 
typhoon  intensity. 

Another  example  of  the  quality  of  the 
NOGAPS  analysis  product  in  the  tropics  is  illus¬ 
trated  in  Figure  10,  which  shows  a  1000  mb  wind 
nowcast  in  the  western  equatorial  Pacific,  where 


FIGURE  10. 1000  mb  wind  analysis  valid  at  0000  UTC, 
March  5. 1991,  showing  the  coupled  tropical  circulations 
that  are  associated  with  the  strong  zonal  westerly  wind 
anomalies  along  the  equator  in  the  western  Pacific. 
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a  tropical  cyclone  couplet  is  developing.  These 
couplets  are  of  interest  because  the  concentra¬ 
tion  of  strong  westerly  winds  along  the  equator 
between  the  two  storms  can  generate  ocean 
KeK'in  waves  that  may  play  a  role  in  the  El  Nino 
warming  in  the  central  and  eastern  equatorial 
Pacific  (Keen,  1982).  The  nowcast  shown  in 
f^igure  10  is  valid  at  0000  UTC,  March  5,  1991.  In 
the  Southern  Hemisphere,  Kelvin  reached  tropi¬ 
cal  storm  strength  on  March  3.  Since  that  time, 
the  NOGAPS  analysis  used  synthetic  soundings 
to  position  and  enhance  Kelvin’s  circulation. 

I  lowever,  at  the  time  of  this  nowcast,  the 
Northern  Hemisphere  storm  was  just  beginning 
to  develop.  It  did  not  become  tropical  storm 
Sharon  until  0600  UTC  on  March  7.  Accordingly, 
synthetic  observations  for  Sharon  were  not  used 
prior  to  that  time.  Even  so,  the  equatorial  west¬ 
erly  wind  anomaly  associated  with  these  two 
stonns  is  well  depicted.  This  westerly  wind 
event  was  verified  by  independent  island  reports 
that,  at  the  time,  were  not  included  in  the 
NOGAPS  analysis  product.  The  ability  of  an 
atmospheric  model  to  depict  such  events  in  their 
fonnative  stages  may  be  extremely  useful  to  cli¬ 
matologists,  given  the  typical  sparseness  of  data 
in  the  region  of  interest  in  the  western  Pacific 
(Harrison  and  Giese,  1991). 

SUMMARY 

he  Navy’s  operational  global-data  assimila¬ 
tion  system  has  shown  a  notable  improve¬ 
ment  over  the  past  few  years,  and  each  compo¬ 
nent  of  the  data  assimilation  system  has  con¬ 
tributed  to  the  overall  performance.  Recent 
improvements  in  the  data  analysis  component 
include  a  better  specification  of  the  observation 
errors  for  marine  surface  data  and  the  inclusion 
of  new  sources  of  data,  including  SSM/I  and 
Pacific  island  winds  and  synthetic  tropical 
soundings.  As  we  continue  to  monitor  the  per¬ 
formance  of  the  multivariate  optimum  interpola¬ 
tion  analysis,  we  expect  to  further  fine-tune  the 
analysis  statistics.  We  are  also  prepared  to 
assimilate  new  data  .sources  as  they  become 
available.  Finally,  improved  coupling  of  the 
oceanographic  and  atmospheric  data  assimila¬ 
tion  systems  will  continue  to  be  at  the  forefront 
of  future  research  efforts. 
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